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Abstract 

A complete theoretical model of the nuclear spin conversion in 13 CH3F in- 
duced by intramolecular ortho-para state mixing is proposed. The model 
contains parameters determined from the level-crossing spectra of the 13 CHaF 
spin conversion. This set of parameters includes the ortho-para decoherence 
rate, the magnitude of the hyperfine spin-spin interaction between the molec- 
ular nuclei and the energy gap between the mixed ortho and para states. 
These parameters are found to be in a good agreement with their theoretical 
estimates. 
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I. INTRODUCTION 



The study of nuclear spin isomers of molecules was started by the discovery of the 
ortho and para hydrogen in the late 1920's Jl[|. It became clear already in that time that 
many other symmetrical molecules should have nuclear spin isomers too. Nevertheless, 
their investigation has been postponed by almost 60 years. The reason for this delay was 
severe difficulties in the enrichment of spin isomers. The situation is improving now (see 
the review in Ref. 0) but yet we are at the very early stage of this research: in addition 
to the well-known spin isomers of H 2 only a few molecules have been investigated so far. 
Among them, the CH 3 F nuclear spin isomers occupy a special place being the most studied 
and understood. 

The conversion of CH 3 F nuclear spin isomers has been explained H|J in the framework 
of quantum relaxation ||. which is based on the intramolecular ortho-para state mixing 
and on the interruption of this mixing by collisions. This mechanism of spin conversion has 
a few striking features. The nuclear spin states of CH 3 F appeared to be extremely stable 
surviving 10 9 — 10 10 collisions. Each of the collision changes the energy of the molecule 
by 10 — 100 cm -1 and shuffles the molecular rotational state substantially. Nevertheless, 
the model predicts that the spin conversion is governed by tiny intramolecular interactions 
having the energy ~ 10~ 6 cm -1 . 

Under these circumstances, the validity of the proposed theoretical model should be 
checked with great care. This is especially important because the CH3F case gives us the 
first evidence of the new mechanism behind the nuclear spin conversion in molecules. Hy- 
drogen spin conversion, which is the only other comprehensively studied case, is due to the 
completely different process based on direct collisional transitions between ortho and para 
states of H 2 . 

Presently there is substantial amount of the experimental data on CH 3 F isomer conver- 
sion (see || and references therein). Theory and experiment on the CH 3 F isomer conversion 
were compared in a number of papers but these comparisons were never aimed to determine 
a complete set of parameters necessary for a quantitative description of the process. The 
purpose of the present paper is to construct such self-consistent theoretical model of the 
CH 3 F isomer conversion. 
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II. QUANTUM RELAXATION 



The CH 3 F molecule is a symmetric top having the C3„ symmetry. The total spin of the 
three hydrogen nuclei in the molecule can be equal to I = 3/2 (ortho isomers), or / = 1/2 
(para isomers). The values of the molecular angular momentum projection on the molecular 
symmetry axis (K) are specific for these spin isomers. Only K divisible by 3 are allowed 
for the ortho isomers. All other K are allowed for the para isomers. Consequently, the 
rotational states of CH3F form two subspaces which are shown in Fig. [l] for the particular 
case of the 13 CHsF molecules. 

Let us briefly recall the physical picture of the CH 3 F spin conversion by quantum re- 
laxation. Suppose that a test molecule was placed initially in the ortho subspace of the 
molecular states (Fig. [I]). Due to collisions in the bulk the test molecule will undergo fast 
rotational relaxation inside the ortho subspace. This running up and down along the ortho 
ladder proceeds until the molecule reaches the ortho state m which is mixed with the para 
state n by the intramolecular perturbation V. Then, during the free flight just after this 
collision, the perturbation V mixes the para state n with the ortho state m. Consequently, 
the next collision is able to move the molecule to other para states and thus to localize it 
inside the para subspace. Such mechanism of spin isomer conversion was proposed in the 
theoretical paper ||. 

The quantum relaxation of spin isomers can be quantitatively described in the framework 
of the kinetic equation for density matrix Let us consider first a free molecule which is 
not subjected to an external field. One needs to split the molecular Hamiltonian into two 
parts 

H = H + hV, (1) 

where the main part of the Hamiltonian, H , has pure ortho and para states as the eigen- 
states; the perturbation V mixes the ortho and para states. In the first order perturbation 
theory the nuclear spin conversion rate, 7, is given by 

9F IT/ I 2 

7= E Isi J (W p (a') + W (a)), (2) 
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where T a i a is the decay rate of the off-diagonal density matrix element p a i a {a 1 G para; a G 
ortho); hu a > a is the energy gap between the states a' and a; W p (a') and W (a) are the 
Boltzmann factors of the corresponding states. The parameters T a > a , V a > a , and oo a ' a are 
crucial for the quantitative theoretical description of the 13 Cfi3F spin isomer conversion. 

All previous comparisons between the experiment on the CH3F spin conversion and the 
theory were performed using "total" rates of conversion which summarize all contributions 
to the rate from many ortho-para level pairs. The "total" rate is just a single number and 
obviously cannot provide unambiguous determination of all parameters which are present in 
the expression (Q). One may combine the experimental data on "total" rates with theoretical 
calculations of some parameters but it is not easy. In this case one has to perform extensive 
calculations of the intramolecular ortho-para state mixing. Even more difficult is to calculate 
the decoherence rates T a / a . Consequently, development of the self-consistent model of the 
nuclear spin conversion in which all parameters are unambiguously determined should be 
based on a different approach. 

III. LEVEL CROSSING RESONANCES 

Theoretical model of spin conversion predicts strong dependence of the conversion rate, 
7, on the level spacing u a i a (see Eq. @). This can be used to single out the contribution 
to the conversion from each level pair which should substantially simplify the quantitative 
comparison between theory and experiment. It was proposed in [7] and performed in || to 
use the Stark effect for crossing the ortho and para states of CH3F. These crossings result in 
sharp increase of the conversion rate 7 giving the conversion spectra if electric field is varied. 
The experimental data [Q are presented in Fig. 0. It is evident that such spectrum contains 
much more information than the "total" conversion rate which is just a single number. 

Comparison of the conversion spectrum in Fig. ||] with the theory needs a modification of 
the model in order to incorporate the Stark effect. Homogeneous electric field lifts partially 
the degeneracy of the a-states of CH 3 F (see Appendix). The new states, /z-basis, can be 
found in a standard way ||: 

\H>= £ > \a F > \a c >; £ = 0,1. (3) 
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Because electric field in the experiment is relatively small, it is sufficient to consider only 
diagonal matrix elements of the Stark perturbation over angular momentum J when calcu- 
lating the //-states. Energy of the //-states are given by the expression 

£(//) = E free (J,K) + i-^^Y^l (4) 

where Ef ree (J, K) is the energy of free molecule; d is the molecular permanent electric dipole 
moment; S is the electric field strength. The new states are still degenerate with respect to 
the spin projections a, op, and crc, and to the sign of M. An account of the Stark effect 
in the spin conversion model is straightforward. Eq. (fj) should be rewritten in the //-basis 
with the level energies determined by the Eq. (H). 



IV. FITTING OF THE EXPERIMENTAL DATA 



Nuclear spin conversion in 13 CH 3 F at zero electric field is governed almost completely 
by mixing of only two level pairs (J'=ll, K'=1)-(J=9, K=3) and (21,l)-(20,3) ||0|. The 
spectrum presented in Fig. ^|is produced by crossings of the M-sublevels of the para (11,1) 
and ortho (9,3) states. This pair of states is mixed by the spin-spin interaction between 
the molecular nuclei Q. There is no contribution to the mixing of this level pair from the 
spin-rotation interaction because of the selection rule for spin-rotation interaction |AJ| < 1 
fTT| . This is fortunate because the spin-spin interaction can be calculated rather accurately. 
Contrary to that, the spin- rotation interaction in CH3F is known only approximately. For 
more details on spin-rotation contribution to the CH 3 F spin conversion see Refs. [JTIl-O] 

The second pair of ortho-para states (21,l)-(20,3), which is also important for the spin 
conversion in 13 CH3F at zero electric field, is mixed by both the spin-spin and spin-rotation 
interactions. The magnitude of the latter is presently unknown. Nevertheless, it does 
not complicate the fitting procedure because in the vicinity of the (11,1)— (9,3) resonances 
presented in Fig. ^|, the (21,l)-(20,3) pair gives very small and almost constant contribution. 

Let us find out now an analytical expression for modelling the experimental data. We 
start by analyzing the contribution to the conversion rate produced by the level pair (11,1)— 
(9,3) which will be denoted as 7 a (£). This contribution can be obtained using the results of 
Refs. 11: 
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7.(£) = E 2 ^f M|2 (W p (aO + WM) ; 

. J' J 2 \ V ^ J 2 . . 

|1/ M , M | 2 = (2J+1)(2J'+1) T 2 . (5) 

-K' K K'-K J \ —M' M M'-M 

Here Vm'm =< /i'jV^I/i > are the matrix elements of the perturbation in which only 
M-indexes were shown explicitly; (:::) stands for the 3j-symbol; T is the magnitude of the 
spin-spin interaction. Note, that the selection rules for the ortho-para state mixing by the 
spin-spin interaction result from Eq. (||): |Aff|; |AJ|; \AM\ < 2. In the fitting procedure T 
will be considered as an adjustable parameter. In Eq. (|5|) we have assumed all T M i M being 
equal: Tm'm = r. This property of V is the consequence of the spherical symmetry of the 
media. The decoherence decay rate V is an another unknown parameter which needs to be 
determined. 

The spacing between the M' and M states in an electric filed follows directly from the 
Eq. ® 

( K'\M'\ K\M\ 

where uq is the gap between the states (J',K') and (J,K) at zero electric field. We have 
considered in Eq. (Q) only pairs of states which have £' = £. They are the only pairs which 
contribute to the spectrum in the electric field range of Fig. The level spacing Uq will 
be considered as an adjustable parameter in the fitting. The dipole moment of 13 CH3F 
in the ground state, which is necessary for the calculation of ujm'm(£), was determined 
very accurately from the laser Stark spectroscopy of 13 CH3F and was found equal d = 



1.8579 ± 0.0006 D 16 



At zero electric field the level pair (21,l)-(20,3) contributes nearly 30% to the total 
conversion rate pi. At electric fields, where 7a(£) has peaks, this contribution is on the 



order of 10~ 2 in comparison with 7 Q (£). The first crossing of the pair (21,l)-(20,3) occurs at 
~ 4000 V/cm thus having its peaks far away from the electric field range of Fig. |2|. In the 
electric field range of Fig. ^| (1-1200 V/cm) the contribution from the pair (21,l)-(20,3) is 
changing by 10% only. Consequently, in the fitting procedure the (21,l)-(20,3) contribution 
is assumed to be constant. This quantity will be denoted as 
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To summarize, the function which will be used to model the experimental data is 

1 (S)= la (S)+ 7b . (7) 

This function contains adjustable parameters T, T, u , and %. 

The result of the least-square fit is shown in Fig. by solid line. The error of the 
individual experimental points in Fig. ^| was estimated as 7%. The values of the parameters 
are given in the Table 1, where one standard deviation of statistical error is indicated. 

Electric field in the Stark cell was determined in experiment H by measuring the voltage 
applied to the electrodes and assuming the distance between them equal to 4.18 mm, which 
is the spacer thickness. It was found out after the experiment || was performed that the 
thickness of the glue used to attach the Stark electrodes was not negligible. The updated 
spacing between the electrodes in the Stark cell is I = 4.22 ± 0.02 mm. Such correction of 
the spacing gives 1% systematic decrease of the experimental electric field values given in 
||. This shift is taken into account in Fig. |2|. 



V. THEORETICAL ESTIMATION OF THE PARAMETERS 

Let us compare the parameters obtained in the previous section with their theoretical 
estimates. We start from the analysis of the level spacing ujq- The best sets of the ground 
state molecular parameters of 13 CH 3 F are given in Ref. [17]. The spacing between the levels 
(11,1) and (9,3) is presented in the Table 1 where the set having most accurate molecular 
parameter A was used. The theoretical value appears to be close to the experimental one 
obtained from the spin conversion spectra. The difference between them is 

LJ (exp) - uj (theor) = 1.0 ± 0.3 MHz, (8) 

which is less than 1% in comparison with ujq itself. 

Next, we calculate the parameter T which characterizes the spin-spin mixing of the level 
pair (11,1)— (9,3) in 13 CH 3 F. The spin-spin interaction between the two magnetic dipoles mi 
and m 2 separated by the distance r has the form ||: 

nv 12 = p 12 m® : T ( 12 ) , 

T// 2) = - 3n^ ; P 12 = mi m2/r 3 I W I i2) , (9) 



where B 1 ) and I*- 2 ) are the spin operators of the particles 1 and 2, respectively; n is the unit 
vector directed along r; % and j are the Cartesian indexes. 

For the spin-spin mixing of the ortho and para states in 13 CH 3 F one has to take into ac- 
count the interaction between the three hydrogen nuclei (Vhh), between the three hydrogen 
and fluorine nuclei (Vhf), and between the three hydrogen and carbon nuclei (Hhc)- Thus 
the total spin-spin interaction responsible for the mixing in 13 CHsF is 

Vss = V HH + V HF + V HC . (10) 

The complete expressions for all components of Vss can be written by using Eq. (|9]) for the 
spin-spin interaction between two particles. For example, for Vhf one has 

Vhf = P H FY,I {n) I F l TnF ; n = 1,2,3. (11) 

n 

Here Phf is the scaling factor analogous to P\ 2 in Eq. n refers to the hydrogen nuclei 
in the molecule. 

T can be calculated in a way similar to that used previously j|. It gives 

|T| 2 = 3\Phh1~2^ | 2 + 2\P H f1~2£\ 2 + ^\ Ph 0^2,2 \ 2 ■ (12) 

Here T 2 1 is the spherical component of the second rank tensor T lq calculated in the molecular 
system of coordinates. The superscripts Iq indicate the interacting particles: 1 refers to the 
hydrogen nucleus H^ 1 ) and q refers to the nucleus of H^ 2 \ or F, or C. 

The calculation of T needs the knowledge of the molecular structure. We used the ground 
state structure of 13 CH 3 F determined in [0,0: r CF = 1.390(1) A, r CH = 1.098(1) A, and 
(3(F — C — H) = 108.7°(2). The numbers in parentheses represent the error bars in units of 
the last digit. By using these parameters one can obtain the value of T which is given in 
the Table 1. The difference between the experimental and theoretical values of T is equal 
to 

%x P - %heor = -5.1 ± 0.5 kHz. (13) 
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VI. DISCUSSION 



Small difference between the experimental and theoretical values of ujq unambigiously 
confirms that the mixed ortho-para level pair (9,3)— (11,1) was determined correctly. From 
the spectroscopical data one can conclude that there are no other ortho-para level pairs 



which can mimic the level spacing ujq = 130.99 MHz. It is also true if one takes into account 
even all ortho-para level pairs ignoring the restrictions imposed by the selection rules for the 
ortho-para state mixing. 

The difference between experimental and theoretical values of the level spacing at zero 
electric field, ujq, is only 1.0 ± 0.3 MHz. The main error in theoretical value of ujq is caused 
by the error in the molecular parameter A . It gives nearly half of the error indicated in 
the Table 1. On the other hand, the J and K dependences of the molecular electric dipole 
moment are too small [16j to affect our determination of the theoretical value of ujq. It is 



possible that the experimental value of uo is affected by the pressure shift, which magnitude 
we presently do not know. Further investigations can precise the frequency gap between the 
states (9,3) and (11,1). 

The difference between the experimental and theoretical values of T is rather small 
(~7%) but well outside the statistical error. This difference may originate from our method 
of calculating T u q^ in which we used the molecular structure (bond lengths and angles) 
averaged over ground state molecular vibration. More correct procedure would be to av- 
erage an exact expression for T over molecular vibration. This requires rather extensive 
calculations. 

There are a few contributions to the systematic error of value of 7j§ . The response 
time of the setup used to measure the concentration of ortho molecules (~ 1 sec) was not 
taken into account in the processing of the experimental data. This gives systematic 
decrease in the value of ^l§y Another few percent of the systematic error may appear due 
to the procedure employed in || to find out the conversion rate inside the Stark cell. This 
procedure relies on the ratio of the Stark cell volume to the volume outside the electric 
field. Taking these circumstances into account we can estimate that up to ~10% difference 
between the experimental and theoretical values of T can be explained by the systematic 
errors. Despite this difference, it is rather safe to conclude that our analysis has proven 
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that the levels (9,3) and (11,1) are indeed mixed by the spin-spin interaction between the 
molecular nuclei. It is impressive that the level-crossing spectrum in the 13 CH 3 F isomer 
conversion has allowed to measure the hyperfine spin-spin coupling with the statistical error 
of 0.5 kHz only. 

Comparison between the measured spectrum and the model supports our choice for the 
T M i M being independent on M and M' . Independence of this parameter on M is the direct 
consequence of the spatial isotropy of the media. The independence on M' — M is more 
intricate. This will be discussed in more detail elsewhere. 

The value of T obtained from the fitting procedure, V = (1.9±0.1) • 10 8 s _1 /Torr, appeared 
to be close to the level population decay rate 1.0 • 10 8 s _1 /Torr measured in Ref. [^] for 
the state (.7=5, K=3) of 13 CH3F. The factor 2 difference is not surprising. T refers to the 
decay rate of the off-diagonal density matrix element p M / M between the states (11,1) and 
(9,3) which should be different from the population decay rate. In addition, the rotational 
quantum numbers in these two cases are different too. 

Column designated as 7(0) in the Table 1 gives the rates at zero electric filed. The 
"theoretical value" is the magnitude of 7(0) given by the solid line in Fig. [| The theoretical 
value coincides well with the experimental one from Ref. [2~I|. Finally we would like to 
mention that our analysis of the spin conversion spectrum has allowed to disentangle for the 
first time the contributions to the conversion rate which arise from the mixing of the two 
level pairs: (9,3)-(ll,l) and (20,3)-(21,l). 



VII. CONCLUSIONS 

We have performed the first quantitative comparison of the level-crossing spectrum of the 
nuclear spin conversion in 13 CH 3 F with the theoretical model. This approach has allowed 
to single out the contribution to the spin conversion caused by the mixing of one particular 
pair of the ortho-para rotational states of the molecule and confirmed unambiguously that 
the mechanism of the intramolecular state mixing is the spin-spin interaction between the 
molecular nuclei. 

All important parameters of the theoretical model which describe the nuclear spin con- 
version in 13 CH3F due to the spin-spin mixing of the ortho-para level pair (9,3)— (11,1) are 
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determined quantitatively. These parameters are the decoherence rate, T, the spin-spin 
mixing strength, T, the level spacing, uj , and the contributions to the conversion rate from 
the two level pairs separately (9,3)-(ll,l) and (20,3)-(21,l). While the decoherence rate V 
is difficult to estimate on the basis of independent information, the experimental values for 
the spin-spin mixing, T, and the level spacing, ujq, are found to be close to their theoretical 
values. These results prove that the nuclear spin conversion in the 13 CH 3 F molecules is 
indeed governed by the quantum relaxation. 
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VIII. APPENDIX 



The CH 3 F quantum states in the ground electronic and vibration state can be classified as 
follows |22]-p4|. CH3F is a rigid symmetric top but it is more transparant to take molecular 



inversion into account and classify the states in D^h symmetry group. First, one has to 
introduce an additional (molecular) system of coordinates which has the orientation defined 
by the numbered hydrogen nuclei and z-axis directed along the molecular symmetry axis. 
Next, one introduces the states 



\/3>= \J,K,M > \I,a,K >; K > 0, (14) 

which are invariant under cyclic permutation of the three hydrogen nuclei: Pu3\/3 >= \(3 >. 
In Eq. (|T4|), \J, K, M > are the familiar rotational states of symmetric top, which are 
characterized by the angular momentum (J), its projection (K) on the z-axis of the molecular 
system of coordinates and the projection (M) on the laboratory quantization axis Z. I and a 
are the total spin of the three hydrogen nuclei and its projection on the Z-axis, respectively. 
The explicit expression for the spin states \I, a, K > is given in WI . 
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Permutation of any two hydrogen nuclei in CH3F inverts z-axis of the molecular system of 
coordinates. Consequently, the action of such operation (P23, for instance) on the molecular 
states reads: -P23I/? >= \P >, where (3 = {J, —K, M, I, a}. Note that the complete set of the 
molecular states comprises both (3 and f3 sets. 

Using the states \f3 > and \(3 > one can construct the states which have the proper 
symmetry with respect to the permutation of any two hydrogen nuclei: 

\/3,k> = [1 + (-1) K P 23 ] |/3 >; k = 0,1. (15) 

The action of the permutation of two hydrogen nuclei on the state \(3, k > is defined by the 
rule: P23\/3, k >= (— l) K \f3, k > and by similar relations for the permutations of the other 
two pairs of hydrogen nuclei. 

In the next step, one has to take into account the symmetric (\s — 1 >) and antisymmet- 
ric (\s = >) inversion states. The action of the permutation of the two hydrogen nuclei 
on these states, for example P 23 , reads 

P 23 | s = > = -\s = >; P 23 \s = 1 > = \s = 1 > . (16) 

Evidently, the cyclic permutation of the three hydrogen nuclei of the molecule does not 
change the inversion states. 

The total spin-rotation states of CH 3 F should be antisymmetric under permutation of 
any two hydrogen nuclei, because protons are fermions. Consequently, the only allowed 
states of CH 3 F are \/3, k = s > \s >. 

Finally, the description of the CH 3 F states should be completed by adding the spin states 
of fluorine and carbon ( 13 C) nuclei, both having spin equal 1/2: 

I at > = \/3, K = s > \s > \cjf > \<Jc >■, (17) 

where of and oc are the Z-projections of the F and 13 C nuclei' spins, respectively. In 
the following, we will denote the states (|17D of a free molecule as a-basis. For the rigid 
symmetric tops, as CH 3 F is, the states \a > are degenerate over the quantum numbers s, 
M, a, crp, a c . 
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Table 1. Experimental and theoretical parameters of the nuclear spin conversion in 
CH 3 F by quantum relaxation. 





ujo/2-k 


r 


r 


7(0) 


lb 




(MHz) 


(kHz) 


(10 8 s-VTorr) 


(10~ 3 s^/Torr) 


(10~ 3 s^/Torr) 


Experiment 


132.06 ± 0.27 


64.1 ± 0.5 


1.9 ± 0.1 


12.2±0.6 (4) 


4.6 ± 0.7 


Theory 


130.99 ±0.15^ 


69.2±0.2( 2 ) 


l.O® 


12.04 ±0.5^ 




Difference 


1.0 ± 0.3 


-5.1 ± 0.5 




0.15 ± 0.8 





^ Calculated using the molecular parameters from Ref. ||17|| , (Table 1, column 2). 
^ Calculated using the molecular structure determined in Ref. f!l|[T9"|j . 
^ 3 ^The level population decay rate from Ref. |20 |. 
^ Experimental value from Ref. [21] 



( 5 ) Zero-field value predicted by the theoretical curve in Fig. 0. 
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FIGURES 




FIG. 1. Ortho and para states of CH3F in the ground vibrational state. The level energies 
were calculated using the molecular parameters from [17]. The (ll,l)-(9,3) pair of states is shown 
to be mixed by the intramolecular perturbation. The bent lines indicate the rotational relaxation 
induced by gas collisions. These collisions do not provide direct ortho-para transitions. 
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FIG. 2. Experimental [8] and theoretical ortho-para conversion spectrum in 13 CH^F. 
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